Local-Field Effects in Silicon Nanoclusters 
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The effect of the local fields on the absorption spectra of silicon nanoclusters (NCs), freestanding 
or embedded in SiC>2, is investigated in the DFT-RPA framework for different size and amorphization 
of the samples. We show that local field effects have a great influence on the optical absorption of 
the NCs. Their effect can be described by two separate contributions, both arising from polarization 
effects at the NC interface. First, local fields produce a reduction of the absorption that is stronger 
in the low energy limit. This contribution is a direct consequence of the screening induced by 
polarization effects on the incoming field. Secondly, local fields cause a blue-shift on the main 
absorption peak that has been explained in terms of perturbation of the absorption resonance 
conditions. Both contributions do not depend neither on the NC diameter nor on its amorphization 
degree, while showing a high sensitivity to the environment enclosing the NCs. 
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I. INTRODUCTION 

The indirect nature of the energy band-gap in silicon 
has always been the major obstacle for its employment in 
light emitting devices, since the momentum conservation 
requires additional mechanisms involved in the recom- 
bination process (e.g. electron-phonon interaction), 
that occur with low probability, hence producing very 
poor emitting rates. In the last decade the discovery of 
efficient visible photoluminescence (PL) and optical gain 
from silicon nanoclusters (Si-NCs) has demonstrated the 
possibility to overcome the indirect band gap of silicon 
by exploiting the quantistic behaviour of the matter at 
the nanoscale. 1,2 

Theoretically, the optical emission has been attributed 
to transitions between states localized inside the 
nanocrystal [as a consequence of the so-called quantum 
confinement (QC) effect], 3-6 or between defect states. 7-12 
While there is still some debate on which of the above 
mechanisms primarily determines the emission energy, 
some recent works have proposed that a concomitance of 
both mechanisms is always present, favouring one or the 
other depending on the structural conditions. 13-22 In the 
attempt of explaining the observations, it was suggested 
that for NC diameters above a certain threshold (of 
about 3 nm) the emission peak simply follows the QC 
model, while interface states would assume a crucial role 
only for small-sized NCs. Anyway, such consclusion was 
still unsatisfactory in many case. More recently, in a 
brilliant experiment Godefroo et al. solved the puzzle 
demonstrating, by using magnetic fields to tune the QC 
and UV lasers to induce defects, that it is possible to 
reversibly control the origin of the PL by introducing or 
removing defects in a single sample: in the former case 
the PL originates from defects while in the latter case it 
originates from QC. 

Previous works already highlighted the dramatic sen- 
sitivity of the opto-electronical properties to the Si/Si02 



interface configuration, especially for very small NCs 
(d < 1 nm), where a large proportion of the atoms is 
localized at the interface. In the latter case, several NC 
characteristics such as passivation, symmetry, and strain, 
considerably concur for the determination of the final 
opto-electronic response, producing sensible deviations 
from the simple QC model. 20 Moreover, many PL ex- 
periments demonstrated that only a very small fraction 
of the NCs in the samples contributes to the observed 
PL, enforcing the idea that precise structural conditions 
are required in order to achieve high emission rates. 24-27 
Finally, recent calculations reported especially high opti- 
cal yields for small NCs, 28 , enhancing the weight held by 
their contribution in real samples. It is therefore clear the 
importance of understanding the factors that, at these 
sizes, contribute to enhance (or reduce) the NC optical 
response. 

Embedding Si-NCs in wide band-gap insulators is 
one way to obtain a strong QC. Si-NCs embedded in 
a silica matrix have been obtained by several tech- 
niques as ion implantation, 9 ' 21 ' 22 ' 29 chemical vapour 
deposition, 19 ' 23 ' 30,31 laser pyrolysis, 8,32 electron beam 
lithography, 33 sputtering, 10 ' 34 and some others. 
Experimentally, several factors contribute to make the 
interpretation of measurements on these systems a dif- 
ficult task. For instance, samples show some dispersion 
in the NC size, that is difficult to be controlled. In this 
case it is possible that the observed quantity does not 
correspond exactly to the mean size but instead to the 
most responsive NCs. 25 Moreover, NCs synthesized by 
different techniques often show different properties in 
size, shape and in the interface structure. Finally, in 
solid nanocrystal arrays some collective effects caused 
by electron, photon, and phonon transfer between the 
NCs can strongly influence the electron dynamics in 
comparison with the case of isolated NCs. 9 In practice, 
all the conditions remarked above lead to measurements 
of collective quantities, making the identification of the 
most active configurations at the experimental level a 
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TABLE I: Structural characteristics for the crystalline embedded NCs (top set) and amorphous embedded NCs (bottom set). 
For each structure are reported, respectively: number of Si atoms forming the NC (NC-Si), number of Si atoms forming the 
NC and not bonded with oxygens (core-Si), whether the NC is centered or not on one silicon ( Si- centered) , number of oxygens 
bonded to the NC (interface-O), number of oxygens bridging two NC-Si (bridge-bonds) , average diameter d, supercell volume 
V a . 



non trivial task. 

From the theoretical side, the possibility of atomically 
manipulating the structures and of associating the 
selected configuration to the calculated response allows 
in principle to elucidate some of the points raised above. 
However, it must be taken into account that an accurate 
characterization of the electronical properties requires 
a full ab-initio approach, limiting the systems size to a 
few thousands of atoms in the case of density functional 
theory (DFT) methods. In addition, the calculation 
of realistic optical absorption or emission spectra, 
involving excited states, requires refined treatments 
that dramatically increase the computational effort, 
furthcrly reducing the maximum manageable system size. 

In this work we present DFT calculations in the local 
density approximation (LDA) of the ground-state elec- 
tronic configuration of Si-NCs of different size embedded 
in a SiC>2 matrix, both in the amorphous and in the 
crystalline phase. The reason for taking into account 
amorphous NCs is based on the fact that real samples 
are always characterized by a certain amount of amor- 
phization, in particular for NCs of small diameter. 16,35 > 36 
The absorption spectra, represented by the imaginary 
part of the dielectric function, are then evaluated within 
the random- phase approximation (RPA), with and 
without the inclusion of local field effects (LFE). 
It is well known that the DFT-LDA severely underesti- 
mates the band gaps for semiconductors and insulators. 
A correction to the fundamental band gap is usually 
obtained by calculating the separate electron and hole 
quasiparticle energies via the GW method. 37 In this 
method the self-energy S is expanded in terms of the sin- 
gle particle Green function G and the screened Coulomb 
interaction W, and at the first order is truncated to the 
first term S ~ iGW . The knowledge of the quasiparticle 
energies, however, is still not sufficient to correctly de- 
scribe a process in which electron-hole pairs are created 
such as the light absorption process. If the electron 
and the hole, created during the absorption process, are 
considered as independent quasiparticles, the structures 
of the absorption spectrum are located at the differences 
between the corresponding singlc-quasiparticles excita- 



tion levels. However, the attractive interaction between 
the positively and negatively charged quasiparticles 
can lead to a strong shift of the peak positions as well 
as to distortions of the spectral lineshape, known as 
excitonic effects. Within the many-body perturbation 
theory (MBPT) framework such interaction is taken into 
account by the solution of the Bcthe-Salpeter equation 
for the polarizability. 38 

An alternative approach to MBPT for the computation 
of neutral excitations is represented by time-dependent 
DFT (TDDFT). 39 TDDFT is expected to be more 
efficient than the MBPT-bascd approach, however 
many conceptual and computational problems remains 
still unsolved preventing its application to complex 
systems. 37 Moreover, a recent comparison of the two 
techniques applied to Si-NCs revealed that TDDFT does 
not take into account correctly the screened Coulomb 
interaction, also for small NCs. 40 

Even if complex MBPT treatments should be invoked 
to include self-energy and excitonic effects, previous 
many-body calculations on Si-NCs reported fundamental 
gaps 41 and absorption spectra 40,43 ' 44 very close to the 
independent-particle calculated ones when LFE are 
neglected. In Ref. 43 we verified this statement for the 
Siio and a-Siio embedded NCs, showing that self-energy 
corrections (calculated through the GW method) and 
electron-hole Coulomb corrections (calculated through 
the Bethe-Salpeter equation) nearly cancel out each 
other (with a total correction to the gap smaller than 
0.2 cV). These considerations justify our choice of 
DFT-LDA for the calculation of the optical properties 
of larger clusters, allowing a good compromise between 
results accuracy and computational effort. 



II. STRUCTURES AND METHOD 

The crystalline embedded structures have been obtained 
from a betacristobalite cubic supercell by removing all 
the oxygens included in a cutoff-sphere, whose radius 
determines the size of the NC. By centering the cutoff- 
sphere on one silicon or in an interstitial position it is 
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possible to obtain structures with different symmetries. 
To guarantee a proper shielding of the strain arising from 
the difference in the silicon/silica lattice spacing, we pre- 
served a separation of about 1 nm between the NCs 
replica. 

The glass models have been generated using classi- 
cal molecular dynamics simulations of quenching from 
a melt, followed by ab-initio relaxations (See Ref. 43 
for further details). The amorphous dot structures have 
been obtained, like in the case of the crystalline systems, 
by applying the cutoff sphere on the glass supercells. 42 
It is worth to stress that also in the crystalline case the 
embedding matrix looses its symmetry after the inclu- 
sion of the NC, due to the metastable nature of the 
betacristobalite. 43 

The relaxation of all the structures have been performed 
using the SIESTA code 45 ' 46 with a DZP basis set (double- 
ts basis plus polarization orbitals) and Troullier-Martins 
pscudopotentials with non-linear core corrections. A cut- 
off of 250 Ry on the density and no additional external 
pressure or stress have been applied. Atomic positions 
and cell parameters have been left totally free to move. 
Following the procedure described above we have built 
three crystalline embedded nanostructures, Siio, S117, 
Si32, and their respective amorphous counterparts, a- 
Siio, a-Sii7, a-Si32- The structural characteristics of all 
the systems are reported in Table I. 
In order to investigate the role of the embedding medium 
on the LFE we have also built the freestanding counter- 
parts of the crystalline and amorphous Si32 NCs by ex- 
tracting the Si-NCs from the silica, including or not the 
first shell of interface oxygens, and then by passivating 
the dangling bonds with hydrogen. In this way we have 
obtained four freestanding NCs: Si32-(OH)56, Si32-Hs6, 
a-Si32-045-H42, a-Si32-H4g. For such structures, in order 
to preserve the strain induced by the silica matrix on the 
NCs, 20 we performed a structural relaxation on the sole 
hydrogens, while holding the position of the other atoms. 
The structural properties of the freestanding structures 
are therefore identical to the embedded ones. 
For each structure we have calculated the eigenval- 
ues and eigenfunctions of the Kohn-Sham Hamilto- 
nian using the ESPRESSO package. 47 Calculations have 
been performed using norm-conserving pscudopoten- 
tials within the Local Density Approximation (LDA) 
with a Ceperley- Alder exchange-correlation potential, as 
parametrized by Perdew-Zungcr. An energy cutoff of 60 
Ry on the plane wave basis set has been considered. 
Once the ground-state geometry has been found, the ab- 
sorption spectra have been computed at the DFT-RPA 
level with and without the inclusion of LFE. The absorp- 
tion spectrum, given by the imaginary part of the macro- 
scopic dielectric function £m(w), is connected to the in- 
verse of the microscopic dielectric function Cqq, (q, to) 
through the so called 'macroscopic average': 48 



When LFE arc neglected at the RPA level, 
em(w) = lim^o £oo(g, w) = 1 - linig^o v(q)P§ , 
where v(q) is the Coulomb interaction and P° is the 
irreducible RPA polarizability. This procedure is in 
fact exact in the case of an homogeneous system for 
which the off-diagonal terms of eQQ,(q,to) are null. On 
the other hand, when local fields (LF) are included 
the quantity €qq (q,ui) must be accessed. Very briefly, 
e" 1 is linked to the reducible polarizability \ by the 
relation e _1 = 1 + v\- At the RPA level we have that 
X = P° + P°vx- Hence by calculating P° = -iG°G° 
with G° single particle Green function, we can obtain \ 
and e _1 . 



III. RESULTS 

It is widely known that LFE assume a crucial role for 
systems characterized by strong charge inhomegeneities. 
Instead, for ordered systems like bulk-silicon and be- 
tacristobalite, LFE tend to vanish out. 49 Besides, the 
same rule applies also for completely amorphized sys- 
tems, like silica-glasses, that at last behave as homoge- 
neous materials. 50 In the case of Si/Si02 heterostruc- 
tures, the inhomogeneity is represented by the interface 
that the NC forms with the surrounding silica, and it 
is therefore important to investigate the role of LF for 
systems with different interface conditions. 

In Figure 1 we report the DFT-RPA absorption spectra 
with and without the LFE contribution, for the crys- 
talline and amorphous Siio, Sii7, and Si32 embedded 
NCs. 51 We note that the microscopic field fluctuations 
produce important screening effects on the spectra of all 
the systems, with a damping of the absorption that is 
more effective at low energies. Instead, at energies larger 
than the optical gap of the embedding medium (^5.5 eV 
for DFT-LDA) the NLF and LF spectra present much 
similar profiles, suggesting that in this energy regime 
the absorption is completely due to pure Si02 states, 
for which the LFE are absent. We also observe a blue- 
shift on the main absorption peaks, in particular for the 
smaller systems in the crystalline phase. The origin of 
such a shift is not clear at this point. We can elaborate 
on this aspect with the aid of a modified Lorcntz oscil- 
lator model, adapted in order to include the interface 
polarization effects (IPE) : 

x + 2^uqx + cjqX = Eo sin(cji) — ax . (2) 

In Eq. 2 the oscillator coordinate, x, corresponds to the 
electronic density displacement, coq is the undamped fre- 
quency, £ is a constant (usually called damping ratio) 
associated to self-relaxation processes, E sin(wt) is the 
external field, and ax is the field produced by the inter- 
face polarization, that contrasts the external field. The 
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FIG. 1: (color online) Absorption spectra for the crystalline (top figures) and amorphous (bottom figures) sets of embedded NCs, 
with (dashed line) and without (solid line) the contribution of local fields. 



solution is 



x(t) = E Q 



(2ujlu £) 2 + (lu% + a - lu 2 ) 2 



sin [cot + <ft\ , 
(3) 

where <f> = arctan [2wcjo£ / {co 2 — ui 2 — a)] is the phase- 
shift that determines the screening. The blue-shift 
appearing in the LFE spectra of Fig. 1 emerges in 
the model as a consequence of the additional term, 
—ax, that changes the resonance frequency from loq 
to uj' = y/ujQ + a. In this picture, damping of the 
absorption and shifting of the resonance peak arise from 
the same physical quantity, and are therefore intimately 
connected. 

It is interesting at this point to examine the response 
of the model in upper and lower outermosts of the 
driving field. At low frequencies, w C wj, the interface 
polarization oscillates in phase with the external field 
due to a small phase-shift <fi. The external field is 
therefore maximally screened in this regime. In prin- 
ciple, the same situation occurs at high frequencies, 
w » wj, where the phase-shift is a decreasing function 
approaching zero. Besides, the function of co multiplying 
the sine in Eq. 3 tends to zero in the limit of high 
frequencies, leading to vanishing LFE. At lo ~ uj' 
the interface polarization oscillates tt/2 out of phase 
from the external field that gets therefore anti-screened, 
leading to the formation of the (shifted) absorption peak. 

It is worth to note that LFE have a strong influence at 
every NC size. This result is in agreement with previ- 
ous works reporting high LFE on very large nanostruc- 
tures, suggesting that they arise almost entirely by clas- 



sical effects. 52,53 In addition, we observe similar trends of 
the spectra of the two largest NCs (S132 and a-Si32). This 
supports the idea that for large NCs the response should 
depend mainly on the IPE and not on the amorphization 
degree, nor on the particular geometry of the Si/Si02 
interface. Finally, differently from Ref. 54, we observe 
similar absorption magnitudes for different NCs size (es- 
pecially in the more realistic case of amorphous NCs), in 
agreement with recent experimental observations. 55 
The comparison of the LF spectra of Figs. lc,lf with 
experimental measurements on Inm-sized Si/Si02 NCs 
show a nice match of the absorption profile and of the 
maximum absorption energy located at about 6 eV. 56 
In order to discuss the role of the embedding medium 
we have considered also the case of freestanding NCs. In 
Ref. 43 we showed that, while the freestanding hydro- 
genated NCs present a larger band-gap and miss oxygen- 
related states at the band-edge, the NC+interface system 
is able (when the strain induced by the embedding ma- 
trix is preserved) to nicely reproduce the characteristics 
of the full NC+Si0 2 system. 

In Figure 2 we report the absorption spectra for the set 
of freestanding Si32 NCs, with and without the inclusion 
of the LFE. First of all we note that the NLF results 
confirm the point remarked above: the NLF spectra of 
the OH-terminated NCs (Figs. 2a and 2c) well match 
those of the corresponding embedded NCs in the 0-6 eV 
range, while both the NLF spectra of the hydrogenated 
NCs (Figs. 2b and 2d) present some modifications and a 
blue-shift of about 0.5 cV on the main peak. 
Also the LF spectra of the NC+interface systems 
match nicely those of the corresponding full systems, 
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evidencing a strong screening in the 0-4 eV range and a 
blue-shift of about 2 eV on the main peak. These spectra 
arc nicely comparable with those of Ref. 56 obtained 
by an effective-medium approach from experimental 
measurements on embedded NCs. It is worth to stress 
that in our samples the crystallinity degree and the 
amorphization degree are intentionally maximized. It is 
therefore reasonable to observe spectra of real lnm-sizcd 
NCs lying in-between those of Figs. 2a and 2c. 
In the case of hydrogenated NCs the LFE become 




Energy (eV) 

FIG. 2: (color online) Absorption spectra for the S132 free- 
standing NCs, with (dashed line) and without (solid line) the 
contribution of local fields: (a) crystalline OH-terminated, (b) 
crystalline H-terminated, (c) amorphous OH-terminated, (d) 
amorphous H-terminated. The spectra are in arbitrary units. 



of dramatic proportions, with severe screening effects 
and a shifting of the main peak of about 5 eV, both 
ascribable (by the model discussed above) to a faster 
response of the iterface polarization. At last it is worth 
noticing that, as already evidenced for the embedded 
case, the introduction of LFE produces spectra that are 
poorly dependent on the amorphization degree. This 
consideration, together with the fact that in real samples 
the smaller NCs tend to remain amorphous, 16 ' 35,36 is 
positive in photovoltaics applications, where an efficient 
harvesting of the sunlight spectrum is required. The 
small variation of the response on the amorphization 
degree is in this case an advantage, because any variation 
could rule the smallest NCs out of the absorption process. 

In the last part of this section we investigate on the dif- 
ferent contributions of Si-NC, interface, and enclosing en- 
vironment, to the absorption spectra. 
In Figure 3 we report the density of states (DOS) pro- 
jected on the NC silicons, NC+interface, and all the 
atoms, of the embedded and freestanding Si32 systems. 
All the projected DOS (PDOS) have been normalized 
following the constraint 

f " PDOS(E)dE=l , (4) 

J —OO 

where PDOS(E) is the DOS projected on the NC sili- 
cons, and Ef is the Fermi energy located half-between 
the energy band gap. 

In the case of embedded NCs (Figs. 3a, 3d) the contribu- 
tion due to the interface-oxygens is very small concerning 
the conduction band while it assumes an important role 
in the valence band, especially for energies below -1 eV 
(-2 eV in the amorphous case). It follows that, at least at 
the NLF level, the absorption below ~3.4 eV (-2.9 eV 
in the amorphous case) is mostly due to states localized 
on the NC, while at higher energies also interface-to-NC 
transitions occur. As expected, the PDOS for the OH- 
terminated NCs (Figs. 3b, 3e) closely reproduce those 
of the corresponding embedded counterparts, confirming 
the idea that the NC+intcrface system is non-interacting 
with the remaining atoms. 43,51 Instead, in the hydro- 
genated case the PDOS of the NCs present a raising in 
the -2 - -5 eV region with respect to the other systems, 
that is responsible for the increased absorption in the 6-8 
eV (see NLF spectra of Figs. 2b, 2d). 



IV. CONCLUSIONS 

Local field effects (LFE) have a great influence on the 
optical absorption of silicon nanoclusters (Si-NCs), and 
must be included for a realistic description of the optical 
response. Their effect arise from the polarizations 
occurring at the NC interface (IPE), and they are 
therefore particularly sensitive to the details of the NC 
termination. 
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FIG. 3: (color online) Projected-DOS for the Si?,2 NCs: (a) cristalline embedded, (b) crystalline OH-terminated, (c) crystalline 
H-terminated, (d) amorphous embedded, (e) amorphous OH-terminated, (f) amorphous H-terminated. The projections on the 
NC silicons, NC+interface, and all the atoms are highlighted in orange (gray), light-orange (light- gray) , and black, respectively. 
The PDOS are normalized following the procedure described in the text. A gaussian broadening of 0.1 eV has been applied to 
the data. 



The IPE depend on the energy of the incoming field. 
At low energies the incoming field and the induced 
polarizations are in-phase, resulting in a screened effec- 
tive field and a subsequent damping of the absorption. 
Contrary, at high energies the polarizations are not able 
to follow the fast oscillations of the incoming field, and 
LFE tend to vanish out. Besides, at energies for which 
the incoming field and the induced polarization are op- 
portunely dephased, the former becomes anti-screened, 
leading to a maximized absorption. Therefore, the final 
(corrected) spectrum appears reduced and blue-shifted, 
by magnitudes that depend on the interface polarizabil- 
ity conditions. In presence of interface oxygens, as for 
freestanding like as for embedded NCs, the LFE produce 
a severe reduction of the absorption up to 4 eV, and 
a blue-shift on the main peak settling around 2 eV. In 
the case of hydrogenated NCs the LFE are of dramatic 
proportions, with a severe damping of the absorption 
in the 0-4 eV energy range, and a blue-shift attaining 
an impressive value of about 5 eV. Interestingly, in all 
the case considered, while the uncorrected spectra show 
an important sensitivity to the structural configuration 
of the NCs (amorphization), LFE tend to smooth out 
such differences in favor of a more consistent description 
of the absorption characteristic. This result suggests 



the possibility of relying on simpler methods for the 
evaluation of the LFE, like the effective medium theories 
as already suggested for other silicon nanostructures. 52 
Such a simplification would be especially conve- 
nient when considering larger NCs, requiring to date 
prohibitive computational efforts with the current model. 
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